7 


MQ-A092  669 


UNCLASSIFIED 


STEVENS  INST  OF  TECH  HOBOKEN  NJ  DAVIDSON  LAB  F/G  13/10 

BLADE  PRESSURE  DISTRIBUTION  FOR  A  MODERATELY  LOADED  PROPELLER. (U) 

SEP  80  S  TSAKONASr  J  P  BRESLIN.  W  R  JACOBS  N00014-77-C-0059 

NL 


i) 26  3  9 


s 


DTIC 

ELECTED 

DEC  8  1980 


This  study  was  sponsored  by  the 

Naval  Sea  Systems  Command 
General  Hydromechanics  Research  Program 
Under  Contract  N000-I4-77-C-0059 
Administered  by 

UWT  Naval  Ship  Research  and  Development  Center 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


(DL  Project  4475/009) 

80  12  04  u70 


ist-'Z-1* 


SECURITY  CLASSIC  AYton  OT  THIS  PAGE  rWh*n  Date  Entered) 


EPORT  DOCUMENTATION  PAGE 


4.  Title  (and  Subtitle)  - - 

1  BLADE  PRESSURE  DISTRIBUTION  FOR  A 
MODERATELY  LOADED  PROPELLER, 


RKAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


TYPE  OF  REPORT  &  PERlOO  COVERED 

7j  FINAL 

6.'f-^«1TFOmfTTrC'OR'G.  REPORT  NUMBER 


7.  AUTMOROI 


S./Tsakonas,  J;P.;Bresiin  a&d  W. R. , Jacobs 


a.  contract  or  grant  numbers; 

-N0OO^l4-77-C-OO59^y 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
_ AREA  A  aORS  UNIT  LUMBERS 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Davidson  Laboratory/  f // 

Stevens  Institute  of  Technology  /./- 

Castle  Point  Station,  Hoboken,  NJ  07030 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

David  W.  Taylor  Naval  Ship  Research  and  / 

Development  Center,  Code  1505  _ 

Bethesda,  MD  20084 


M.  MONITORING  AGENCY  NAME  »  ADDRESSf/I  dlltefanl  trom  Controlling  Ottlca)  IS.  SECURITY  CLASS,  (ol  lhl«  r«poftJ 

°f.Nava’  Resfrch  UNCLASSIFIED 

800  N.  Quincy  Street 


0>i~6TT53NLBO23OJ  ; 
’  ySR  Jp23  £1  ,01 


12.  REPORT  DATE 

Sep  tutor  1980 


13.  NUMBER  OF  PAGES 

x  +  59  pp. 


Arlington,  VA  22217 


1  ba.  DECL  ASS|F|CATIOn/  OOWN  GRADING 
SCHEOULE 


16.  DISTRIBUTION  STATEMENT  t ol  thla  Report) 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  antarad  In  Block  20,  ft  dlffarant  from  Report ) 


is.  supplementary  notes 

Sponsored  by  the  Naval  Sea  Systems  Command  General  Hydromechanics  Research 
Program  administered  by  the  David  W.  Taylor  Naval  Ship  Research  and  Develop¬ 
ment  Center,  Code  1505,  Bethesda,  MD  20084. 


19.  KEY  WOROS  (Contfnua  on  revaraa  aide  if  necaaaary  and  IdantUy  by  block  numbar) 

Hydrodynamics 

Moderately  Loaded  Propellers 
Blade  pressure  Distribution 


20.  ABSTRACT  (Continua  on  ravaraa  aide  It  nar.aaaary  and  IdantUy  by  block  numbar) 

^Aa  new  theoretical  procedure  has  been  developed  modifying  the  existing 
analysis  for  a  marine  propeller  operating  in  a  nonuniform  inflow  field  by 
considering  the  radially  varying  mean  wake  and  mean  propeller  induction. 

In  addition,  the  selection  of  a  new  reference  surface  around  which  the 
perturbation  analysis  is  developed  Is  based  on  the  nonlinear  form  of  the 

Bernoulli  equation  together  with  an  appropriate  kinematic  condition  ■— - - 

(Continued) 


DO  1473  EDITION  OF  I  NOV  «  IS  OBSOLETE  IIMPI  A«  IFIFD  1 3 

S/N  0  102-0  14- *601  v  - UNILMSS  1  P  I  - 

y  s  (  j  t  i  <  )  SECURITY  CLASSIFICATION  OF  THIS  PAGE  f#>»en  Data  Bntarad) 

1  r  J  /  /  ' 


_ UNCLASSIFIED _ 

II  V  CUtMUHIInw  Ini'.  >•«<,.  LI  I1.lt  Cnlorcl) 

20.  Abstract  (Cont'd) 

- > 

existing  at  the  propeller  operational  condition.  A  flow  field  closer 
to  the  propeller  operating  condition  is  achieved  thereby  and  the 
linear  theory  requirement  of  small  perturbation  quantities  is  rein¬ 
forced.  The  approach  is  thus  applicable  to  moderately  to  heavily 
loaded  propellers  immersed  in  the  stronger  wakes  of  hulls  of  large 
block  coefficient,  although  it  can  be  used  for  lightly  loaded  pro¬ 
pellers  as  wel 1 . 

7; 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  This  FAeSfHTi.n  £>•!«  Cnl.r.F; 


STEVENS  INSTITUTE  OF  TECHNOLOGY 

DAVIDSON  LABORATORY 
CASTLE  POINT  STATION 
HOBOKEN.  NEW  JERSEY 

Report  SIT-DL-80-9-2063 
September  1 980 


BLADE  PRESSURE  DISTRIBUTION 
FOR  A  MODERATELY  LOADED  PROPELLER 


by 

S.  Tsakonas,  J.P.  Breslin  and  W.R.  Jacobs 


I 

1 


I 


This  study  was  sponsored  by  the 

Naval  Sea  Systems  Command 
General  Hydromechanics  Research  Program 
Under  Contract  N000-14-77"C-O059 
Administered  by 

OWT  Naval  Ship  Research  and  Development  Center 
(DL  Pro^-ct  4475/009) 


I 


Approved; 

John  P.  Breslin 

x  +  59  pp.  Di rector 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


I 


ABSTRACT 


A  new  theoretical  procedure  has  been  developed  modifying  the  exist¬ 
ing  analysis  for  a  marine  propeller  operating  in  a  nonuniform  inflow  field 
by  considering  the  radially  varying  mean  wake  and  mean  propeller  induction. 
In  addition,  the  selection  of  a  new  reference  surface  around  which  the 
perturbation  analysis  is  developed  is  based  on  the  nonlinear  form  of  the 
Bernoulli  equation  together  with  an  appropriate  kinematic  condition  exist¬ 
ing  at  the  propeller  operational  condition.  A  flow  field  closer  to  the 
propeller  operating  condition  is  achieved  thereby  and  the  linear  theory 
requirement  of  small  perturbation  quantities  is  reinforced.  The  approach 
is  thus  applicable  to  moderately  to  heavily  loaded  propellers  immersed  in 
the  stronger  wakes  of  hulls  of  large  block  coefficient,  although  it  can  be 
used  for  lightly  loaded  propellers  as  well. 


KEYWORDS 

Hydrodynamics 

Moderately  Loaded  Propellers 
Blade  Pressure  Distribution 
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INTRODUCTION 

In  the  course  of  a  series  of  investigations  concerned  with  ad- 

i  2* 

aptation  of  the  unsteady  lifting  surface  theory  to  marine  propellers,  ’ 
the  following  baste  assumptions  have  been  made: 

a)  the  propeller  operates  in  a  spatially  nonuniform, 
nonseparating  and  noncavitating  flow  of  an  inviscid 
and  incompressible  fluid, 

b)  all  perturbation  quantities  are  considered  to  be  small, 
so  that  the  linearized  version  of  the  theory  may  be  used 
as  a  basis, 

c)  a  helicoidal  surface  of  constant  pitch  (from  hub-to-tip) 
is  selected  as  reference  surface  and  deviations  from  this 
surface  are  considered  as  the  perturbation  quantities. 

The  reference  surface,  over  which  the  singularities  are  distributed, 
is  the  "zero  lift"  helicoidal  with  constant  pitch  2rr  ^  where  U  is 
the  forward  speed  (ship  speed)  and  Q  the  angular  velocity  of  the 
rotating  propeller.  Any  deviation  from  this  surface  is  taken  to  be 
induced  by  perturbation  velocities  such  as  those  due  to  ship  wake, 
incident  flow  angle,  blade  camber,  non-planar  blade  thickness  and 
flow-distortion  blade  thickness.  These  effects  have  been  considered 
separately  and  then,  as  permitted  in  the  linear  theory,  have  been 
added  together  to  determine  their  combined  effects  on  the  loading  and 
hydrodynamic  forces.  The  linearized  version  of  the  unsteady  lifting 
surface  theory  leads  to  an  integral  equation  relating  the  unknown 
loading  distribution  with  the  known  onset  velocities.  The  equation 
is  valid  for  a  lightly  loaded  propeller  with  wakes  of  small  magnitude 
(intensity).  With  today's  tendency  toward  longer  and  fuller  hull  forms,  the 
wake  intensity  behind  the  hull  is  quite  strong.  Thus  the  perturbation  quan¬ 
tities  are  larger  so  that  the  linearized  theory  is  inadequate  for  the  needs 
of  the  propeller  designer  in  regard  to  the  prediction  of  mean  thrust  and 
torque. 


Superior  numbers  in  text  matter  refer  to  similarly  numbered  references 
listed  at  the  end  of  this  report. 
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The  propeller  operating  in  ship  wake  which  varies  with  radial 
and  angular  position  not  only  encounters  the  spatial  inhomogeneities 
of  the  incoming  flow  but  it  is  also  under  the  influence  of  the  self- 
induced  velocity  field.  The  latter  can  be  determined  provided  the 
blade  loadings  are  known  which,  in  effect,  requires  "apriori"  knowledge 
of  the  reference  surface.  Thus  an  iterative  procedure  must  be  estab¬ 
lished  based,  as  a  first  step,  on  an  assumed  helicoidal  reference 
surface  determined  by  the  known  variable  hydrodynamic  pitch  angle  from 
which  the  first  approx imat i ons  of  blade  loading  and  propeller  induction 
are  determined.  Having  this  induction  and  the  radial  variation  of  the 
speed  of  advance  a  new  reference  surface  is  establ ished  and  thus  another 
set  of  calculations  is  performed  for  the  evaluation  of  the  new  blade 
loading  and  induction  effect.  The  iterative  procedure  is  continued 
until  no  variation  in  the  pitch  of  the  reference  Surface  and  in  the 
blade  loading  can  be  observed. 

While  the  iterative  procedure  uses  the  same  linearized  approach 
as  developed  in  Reference  1,  it  takes  into  account  the  quadratic  form 
of  the  Bernoulli  equation,  since  the  pertubation  velocities  are  some¬ 
times  of  the  same  order  of  magnitude  as  the  velocities  of  the  undisturbed 
flow  in  which  case  linearization  of  the  pressure  equation  is  not  valid. 
This  procedure  is  appropriate  to  moderately  loaded  propeller  opera¬ 
tional  conditions.  Once  the  final  geometry  of  the  reference  surface 
is  established,  the  solution  of  the  integral  equation  relating  the 
unknown  loading  with  the  known  onset  velocities  will  determine  the 
blade  pressure  distribution  and  the  corresponding  hydrodynamic  forces. 

It  should  be  noted  that  the  imposed  boundary  conditions  on  the  blade 
are  fixed  and  do  not  change  during  the  iterative  process. 

The  same  reference  surface  is  further  used  for  the  unsteady  flow 
conditions,  since  the  unsteady  onset  velocities  being  small  compared 
to  steady  state  velocities  can  be  considered  as  perturbation  without 
running  the  risk  of  violating  the  basic  requirements  of  the  linearized 
theory. 
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ANALYSIS 


1)  Integral  Equation 

A  theory  and  computer  program  have  been  developed  for  the  case 
of  a  lightly  loaded  propeller  utilizing  the  small  perturbation  approxi¬ 
mation  by  selecting  as  a  reference  surface  the  helicoidal  surface  of 
"zero  lift"  with  constant  pitch  (i.e.,  P  =  2r\  £  =  —  where  U  =  ship 
speed  and  ft  =  propeller  rotational  velocity).  However,  many  propellers 
operate  in  strong  wakes  and  at  lower  relative  inflow  velocities  resulting 
in  larger  deviation  from  the  reference  surface  hitherto  employed.  These  large 
excursions  violate  the  requirement  of  the  small  perturbation  theory. 

A  new  analysis  is  developed  herein  along  the  same  general 
lines  as  the  small  perturbation  theory  but  with  emphasis  on  the 
selection  of  a  more  appropriate  reference  surface  which  simulates  closely 
the  actual  propeller  operational  conditions. 

The  linearized  unsteady  lifting  surface  theory  for  a  marine  pro¬ 
peller,  with  its  blades  lying  on  a  helicoidal  surface  and  operating  in 
the  nonuniform  flow  of  an  incompressible  ideal  fluid,  is  formulated  by 
means  of  acceleration  potential  method.  It  is  based  not  only  on  a  small 
perturbation  approximation  but,  also,  on  the  assumptions  that  the  pro¬ 
peller  blades  are  thin  and  operate  without  cavitation  and  flow  separation. 

It  may  now  be  considered,  more  appropriately,  that  this  propeller 
operates  in  a  spatially  varying  flow  generated  by  a  strong  hull  wake 
and  also  under  the  influence  of  self-induction.  Both  these  factors 
will  influence  the  selection  of  the  proper  reference  surface  in  the 
steady  state  flow  conditions.  It  is  assumed  that  the  flow  field  in 
which  the  propeller  is  immersed  has  velocity  components  given  in  cylindri¬ 
cal  coordinates  by 

UA(r)  +  u.(x,r);  v.(x,r);  w.(x,r) 


1+ 


(1) 
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where  U,(r)  is  the  local  speed  of  advance  and  u.,  v.,  w.  are  the 
A  ill 

axial,  tangential  and  radial  components  of  the  perturbation  velocities. 

The  relation  between  the  velocity  potential  function  $  and 
the  acceleration  potential  function  ijf  has  been  established  by  solving 
the  linearized  Euler  equation  of  motion  to  obtain 

l  x 

$(x,  r,cp;t)  =  — —  j  if(T,  r,cp;t  -  i-l-I)  dr  (2) 

U(r)  I,  U(r) 

where  U(r)  =  L^o,  r)  =  Ufl(r)  +  u.(o,  r)  the  axial  component  of  the 
resulting  velocity  at  the  propeller  plane. 

It  is  further  known  that  the  pressure  field  generated  by  a 
lifting  surface  S  is  given  by  distributed  doublets  with  axis  parallel 
to  the  local  normal  and  with  strength  equal  to  the  pressure  jump  across 
the  surface.  The  reference  surface  over  which  the  doublets  are  dis¬ 
tributed  is  a  helicoidal  surface  with  variable  pitch  along  the  radius. 
These  facts  considered  with  Equation  (2)  lead  (following  Reference  l) 

to  an  integral  equation  which  relates  the  known  onset  velocities 
with  the  unknown  blade  loading  at  each  frequency  q: 

V(r)e  °  =  JT  Aplq\?,p,9o)  K(r,cpo;p,eQ;q)dS  (3) 

S 

where  the  kernel  function  is  given  by 


K(  r,cpo;p,0o;q) 


1 


N 

■  1  im  £ 
4np  ,U(  r)  6_,°  n=  1 


iq9  x 

n  _  r 

Bn'  J 


iqa(r1(T-x) 


3 

3n 


(I)  -T 


w 


Here 


Or  +  v.( r) 

i 

UA(  r)  +  u . (  r) 
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Q  =  angular  propeller  velocity 

U.( r)  =  speed  of  advance  at  the  propeller  plane  at 
radial  location  r 

6-*o  means  x  -»  cp(y/a(r)  and  5  "*  6 Q/a(p) 

0  =  2tt( n- 1 )/N  n  =  1,2. ..N 

N  =  number  of  blades 

R  =  Descartes  distance  between  the  control  point  (x,  r,cpQ) 

1 

and  loading  point  (5,P,6q)  =  [_(  T-§)£+r2+p2-2rpcos(  ®0-cP0+®n“a(  r)(T-x)J 


and 


_  P 

3n  / - - - - 

J 1  +  a2( p)p2 


S_ 

dn'“ 


_ r _ 

J 1  +  a2( r)r2 


p2  ^ 


(  ,  N  a 
^a(r)i^ 


are  the  normal  derivatives  at  the  loading  and  control  point,  respectively. 
It  should  be  noted  that  the  time  factor  has  been  eliminated  from  both 
sides  of  Eq.  (3)  and  it  is  understood  that  only  the  rea  1  pa  rt  of 
the  solution  is  retained.  A  limiting  process  is  introduced  in  the 
kernel  function  to  avoid  the  mathematical  difficulty  due  to  the  presence 
of  a  high  order  singularity.  The  kernel  function  being  one  of  the  most 
complex  in  the  lifting  surface  theory,  attention  is  given  to  the  numeri¬ 
cal  solution  of  the  integral  equation  by  means  of  a  high-speed  digital 
computer.  The  analysis,  however,  has  been  carried  out  to  the  stage  where 
laborious  computations  can  be  efficiently  performed  by  the  numerical 
procedure. 

It  has  been  assumed  that  the  shape  of  the  chordwi se  loading  is 
the  same  as  has  been  shown  by  Landahl  to  be  the  appropriate  one  in 
a  two-dimensional  flow.  Furthermore,  a  method  called  the  "generalized 
lift  operator"  technique^  is  applied  to  both  sides  of  the  integral 
equation  to  reduce  the  surface  Integral  equation  to  a  line  integral 
equation  along  the  propeller  radius.  Then  by  the  collocation  method 
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the  line  integral  equation  is  reduced  to  a  set  of  algebraic  equations 
in  the  unknown  spanwise  loading  distribution. 

The  following  substitutions  are  made  in  Equations  3  and  k: 

*-^(p>ea)  =  (?,P,eo)  *  pe£  ib/ft 

e  =  <tp  -  eq  cose  o  £  e  ^  n 

o  b  a  a 

r  r 

cp  =  a  -  0,  coscp  0  ^  cp  ^  rr 

o  b  Ta 

and  the  expansion  of  the  inverse  Descartes  distance  is 

*  im[0  -tp  +5  -a(r)(r-x)]  «  i(T-§)k 

—  =  -1  £  e  °°n  f  |  ( |l|p)K  ( lklr)e  dk 

R  tt  J  r/  m' 

m« - ® 

for  p  <  r  (otherwise  p  and  r  in  the  modified  Bessel  functions  must 
be  interchanged).  Here  a  is  the  propeller  skewness  (the  angular  posi¬ 
tion  of  the  midchord  line  from  the  generator  line  through  the  hub  in 
the  projected  propeller  plane),  0^  Is  the  subtended  angle  of  the  pro¬ 
jected  blade  semichord,  0  and  cp  are  anqular  chordwise  locations  of 

a  a 

the  loading  and  control  points,  respectively,  and  the  superscripts  p 

and  r  refer  to  the  values  at  the  loading  and  control  points,  respectively. 

After  the  chordwise  integration  and  application  of  the  lift- 
operator,  the  integral  Eq.  (3)  reduces  to: 


e~iq°  q©T)  =  £  J  L^q,n\p)K^m,n\ r,p;q)dp  (5) 

UA(r)  b  n  =1  p 
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and 


t(m,n)  _ 


K 


■  iqACT 


1mp,U.(r)U(r)  a(p)[l+a2(r)i^  m=-“ 

'  A  m=q+2N 


iiN(l-  44 )aP 

•  {  e  ^  a2(r)£N+^\a(p)a(r)iN+  ^]l(m)(q9j)A(n)  Jq+(l- 


•  lm(a(r)UNIp)  Km(a(r)UNIr) 


.  -iiNAo 

-  —  e 

rr 


,klaTrT  "  aTpT;i  ( lklp)K  (  Iklr) 
'  v  7  m _  m 

k-a(  r)j£N 


,m  ((m"  ^T7T^eb^  An((m'  I(77)(0  dk)‘ 


and  Act  =  ar  -  a*3  . 


Here 


rr 


lm(y  )  =  ^  J  *(m) 


iy  coscp 


dcp 


and 


An(z)  =  i  Js<n)  e 


-  izcosB 


“  sine  de 

a  a 


(6) 


the  details  of  which  are  presented  in  Appendix  A. 

The  integral  equation  (5)  is  solved  numerically  by  the  usual 
collocation  method  with  the  loading  L^,n^(p)  assumed  to  be  constant 
over  each  small  radial  strip.  Then  only  the  kernel  needs  to  be  inte¬ 
grated  over  the  radial  strip. 
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Reference  I  gives  details  of  the  analytical  development  and 
the  different  numerical  procedures  used  to  obtain  the  finite  contri¬ 
butions  of  the  Cauchy-type  singularity  of  the  k-integral  at  k  =  a( r) . 
(m  -  q)  and  of  the  higher-order  Hadamard-type  singularity  when  p  =  r. 

2)  Chordwise  Modes 

The  proper  selection  of  chordwise  modes  is  dictated  by  the 
loading  distribution  on  a  foil  in  two-dimensional  flow  conditions. 

The  same  chordwise  modes  have  been  selected  as  those  of  Reference  2  . 

a)  Unsteady  flow  conditions 

The  unknown  chordwise  modes  are  approximated  by  the  known 
Birnbaum  distribution  which  has  the  proper  leading  edge  singularity 
and  satisfies  the  Kutta  condition  at  the  trailing  edge. 

nmax 

L(q)(p,ea)=  .2  L(q'n)(p)  e(n) 

n=l 


n  max 

=  ~  jl(q' 1 )(p)cot  ^  L(q'n)(P)sin(n-1)ea 

(7) 

b)  Steady  state  flow  conditions 

The  unknown  chordwise  mode  shapes  are  selected  to  conform 
to  the  observed  pressure  distribution  on  the  NACA  foil  sections. 

i)  The  NACA-a  mean  line  at  the  design  condition 


The  section  has  a  constant  loading  distribution  from 
—  =  o  at  the  leading  edge  to  —  =  a  (a  varying  from  0  to  I;  x  =  location 
of  points  along  the  chord  and  c  =  chord  length)  and  then  decreases 
linearly  to  zero  at  ^  =  l,  the  trailing  edge.  This  type  is  designated 
as  "roof-top"  loading  and  the  distribution  is  given  by 
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with  the  corresponding  '(z)  given  in  Appendix  A. 

ii)  For  any  other  sections  such  as  NACA  mean  lines  of  the 
k-  and  5-digit  wing  series,  including  the  lenticular  mean  line,  and 
arbitrary  mean  lines  in  general  the  distribution  is 

nmax 

°\  0  )  =  “f  _£  L^°,n\p)  0(n) 

n=l 

nmax 

=  l  jL(o,l)(p)cot  %  +  L(°'n)(p)sin(n-l)eaJ  (9) 

iii)  Off-design  flow  condition 

If  the  solution  of  the  integral  equation  at  the  design 
flow  condition  is  known,  (j^  =  advance  ratio  at  the  design  condition) 
then  the  off-design  condition  can  be  obtained  through  the  input  arrange¬ 
ment.  In  the  steady  state  flow  conditions  at  off-design  advance  ratio 
Jq{j,  the  propeller  is  subjected  to  a  change  in  angle  of  attack  due  to 
AJ  =  JQ(j  -  J^.  The  additional  loading  due  to  this  additional  change 
of  angle  of  attack  will  be  obtained  by  utilizing  as  chordwise  loading 
distribution  the  first  term  of  the  Birnbaum  distribution,  i.e., 

l(o)(P'V  ^l(o)  (p)  cot-r  o°) 

If,  however,  the  information  on  the  design  flow  condition  is 
not  available  then  the  general  form  of  the  Birnbaum  distribution 
should  be  used  for  the  chordwise  loading  distribution 

l(o)(P'V  =  ^  1l(0,,)(p)  cot  T  +  *  L(o,")(p)sin(;-i)e  f 

'  n=  1  ’ 

00 
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3)  Perturbation  Velocity  Distributions 

The  left-hand  side  of  the  integral  equation  represents  the  normal 
component  of  the  velocity  perturbations  from  the  basic  flow  due  to  non¬ 
uniformity  of  the  flow  (wake),  to  blade  camber,  non-planar  blade  thickness 
and  the  incident  flow  angle  in  the  design  cond i t ions,  and  to  an  addi¬ 
tional  flow  angle  in  the  off-design  condition  arising  from  the  difference 
between  off-design  and  design  advance  ratios. 

a)  Normal  velocity  due  to  wake 

The  normal  component  of  the  wake  velocity  along  the  middle 
chord  of  each  radial  strip  is  given  by(cf  Reference  2  ) 

For  q  =  0 


VN(0  V  (r)  V  (r) 

UJTT  =  uTTT  cosep(r)  "  uJTT  sin0p(r) 


where 


Vjr)  Vx(r) 


Vr) 


=  measured  axial  velocity  normalized  with  respect 
to  the  speed  of  advance 


VT(r) 


■jj-»  =  measured  tangential  velocity  normalized  with  respect 

UA'  '  to  the  speed  of  advance 

0  =  blade  pitch  angle  =  tan”'  — 1— ^ 

p  r  a  2nr 

For  a  single  screw  ship  V^.°^  =  0  and  taking  \/r^  =  U.(  r)  then 

l  A  n 


Vr) 

Thus  there  is  no  wake  distribution  in  the  steady  state  case  when  the 
reference  surface  is  selected  on  the  basis  of  the  velocity  of  advance. 
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For  q  ^  0 

The  input  to  the  program  is 

v(q) 

(-TFt)  =  aq(r)  cosep(r)  -  Aq(r)  Sin9p(r) 

UA  R 

V(q\ 

(-Jtj)  =  bq(r)  cosep(r)  -  Bq(  r)  sin0p(r)  (12) 

UA  I 

where  R  and  I  represent  the  real  and  imaginary  parts.  The  coeffi¬ 
cients  a,  b,  A  and  B  are  obtained  from  the  harmonic  analysis  of  the 
wake  survey; 


(r) 

X71 


£  [a  (  r)  cosq® 
q=l  q 


bq(  r> 


sinq©  ] 


V  <r)  . 

— Ir—x  =  X  [a  (r)  cosq® 

V '  q 


+  Bq( r)  sinq®  ] 


The  same  sign  convention  as  in  Reference  2  has  been  utilized.  The 
positive  ®  is  defined  in  the  counterclockwise  direction  from  ©  =  0 
at  the  upright  position  of  the  blade.  The  axial  component  Vx  is 
positive  downstream  and  tangential  component  V^.  is  positive  in  the 
counterclockwise  direction  looking  upstream. 

b)  Normal  velocity  due  to  blade  camber 

The  velocity  V c  normal  to  the  blade  in  the  negative  direction 
induced  by  the  flow  disturbance  caused  by  blade  camber  is  given  by 


V(c°}(r)  ^+a|(r)rfi  <*c(r,s) 

UA(  r)  C(r)  3s 

2-/l+a|(r)r2  3fc(  r,cpQi) 
C(r)sincpa 


(13) 
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where 


Sfc( r,s) 
3s 


slope  of  the  camber  line  f^.(r,s)  measured  from  the 
face  pi tch  1 i ne 


s  =  !( l-cos9&)  chordwise  location  non-dimens ional  ized  on  the 
bas  is  of  C(r) 

C( r)  =  chord  i n  feet 

aft(r)  =  nr/UA(r) 


On  applying  the  lift-operator  to  the  above  expression,  then 


v£0,(r)  2yita=(r)r2  "  9fc  <% 

yTJ  «(T)  {  SM  vr.  ns?, 

where  §  (m)  =  lift  operator  defined  in  Appendix  A. 


(1*0 


The  integration  is  done  as  in  Reference  4.  In  the  small  range 
near  the  leading  edge  0  £  cp^^cos  '(0-9)  or  0  s  s  £  0.05  the  camber¬ 
line  is  assumed  to  be  parabolic  and  in  the  range  cos”  (-.8)  £  cp^  =»  tt 
or  0.90  ^  s  £  1.0  near  the  trailing  edge,  the  camberline  is  assumed 
to  be  straight  line.  The  integration  is  done  analytically  in  these 
regions  and  numerically  over  the  remainder  of  the  chord. 


c)  Normal  velocity  due  to  "non-planar"  blade  thickness 


The  thickness  of  propeller  blades  describing  a  helicoidal 
surface  or  non-planar  surface  generates  a  velocity  field  on  the  blade 
itself  and  consequently,  affects  the  loading  distribution  in  the 
steady  state. 


Resorting  to  the  "thin  body"  approximation,  the  velocity  potential 
due  to  the  blade  thickness  of  an  N-blade  propeller  is  given  by 


*T(x,r,<P0)  =  - 


N 


T-  Z 
^  n=l  § 


p  r  K?,p*e0)  «/i+aA(p)p 

J  l  — 


aA(p)p 


pdpde 
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where 


— - (P,s) 

M(^,P,eo)  =2Ua(p)  Vl+aA(p)p2  - ^ 

sfT(?,P,en) 

,2Uft(P)  — i—  — 


fT(?,p,6Q)  =  thickness  distribution  over  one  side  of  the  blade 

section  at  radial  position  p  in  the  propeller  plane. 

R  =  C(x-§)2  +  r2  +  p2  -  2rp  cos(  e0_cP0+®n)3 2 

e  =  EULlllLl  n  =  1,2  ...  N 
n  N 

0Q,cpo  =  angular  position  of  the  loading  and  control  points, 
respectively,  from  the  reference  line  (see  Figure  1) 

The  corresponding  normal  velocity  component  is  given  by 


v,(o)(0 

UA<r> 


UA( r)  dn‘ 


The  application  of  the  lift  operator  after  the  use  of  the 


transformat  ions 


x  =  — r— \  =  (a  -  e[  coscp  )/aft(  r)  0  £  cp  £  tt 

aA(r)  v  b  Ta''  A'  '  Ta 

e 

F  =  — 7^T  =  (ctP  -  0P  COS0  )/a  (p)  0  s  0  s  it 

b  aA(p;  '  b  or'  Avr/  a 

and  the  expansion  of  1/R  in  a  series  yields  (see  Appendix  B) 

V^°^(r)  - 

*T  '  ’  _  Nr  C  ..  /,  ,  _2,_n  .2 


TT2UA(r)»/l+a2(r)rs  P 


J  UA(p)yi+as(p)p2 


r®  "  fml,  r  o  i u (or/a  ( r) -oP/ a  ( p) 

UaA(r)ulo(up)Ko(ur)L,m'Partll  \"ueb/aA^rVAl\ueb/aA^PVe 


J  1  .  (up)Kje  (ur->  [G(u,£)  +G(u,-X)]du 
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where 


G  ( u  *  -0  =[_aA  (  r)  u+  ^|]lm.Part[l(m)((aA(r)XN-u)e^/aA(r)jA((u-aA(p)£N)eP/aA(p)) 

-UN(crr-afi)  Uj(ar/aA( r)'  °P/aA( P»  -i 


/  \  V  «XCOS0a 

A  X  =  |  e  d0c 


which  can  be  evaluated  once  the  thickness  distribution  f(T)  is  supplied. 

For  any  propeller  the  blade  semi-thickness  distribution  in  the 
projected  plane  can  be  approximated  by  (see  Reference  2) 

fT(P»ea)  -  2p0b  V2P0(p)  sin  +  "TCpT  La0(p)  +  I]  an(p)cos  }  (,8' 

where  fT  and  p  are  fractions  of  propeller  radius  rQ 

PQ(p)  =  ratio  of  leading  edge  radius  to  chord  length  C(p) 
at  given  radial  position 

t^C  =  ratio  of  maximum  thickness  to  chord  at  that  radial 
pos i tion 

The  coefficients  aQ(p)  and  an(p)  are  determined  from  the  equivalent 
fourth  order  polynomial  in  x  by  the  "least  squares"  method  as  shown 
in  Reference  2  . 

d)  Normal  velocity  due  to  flow  angle. 

The  incident  flow  angle  has  been  defined  as  the  difference 
between  the  geometric  pitch  angle  0  ( r)  =  tan  '  and  the  hydro- 

dynamic  pitch  angle  p(  r)  =  tan-^  | yp  of  the  reference  helicoidal 

surface.  The  normal  velocity  due  to  this  effect  is  then  given  by 
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,(°)f 

'f 


vrV) 


UA(r) 


=  “  J  1  +  a2  (r)r2  ^  ep(r  )  -  tan"  1 


which  after  the  application  of  the  lift  operator  becomes 

-1  1 


r,(  °>  m) 
Vf 


-  J\  +  a^(  r)  r2  (^6^  -  tan 


aA(>)r 


for  m=  1  or  2 


for  m  >  2 


(  19) 


since  l^\o)  =  1^2^(0)  =  1  and  1^^(0)  =  0  for  m  >  2. 


At  off -design  advance  ratio  Jod'  the  propeller  is  subjected  to 
a  change  in  angle  of  attack  due  to  AJ  =  JQd  -  J^.  The  normal  velocity 
perturbation  due  to  this  change  in  angle  of  attack  is  given  by 


AV 


rX°,™) 


Uod<r) 


=  {  -  Vi  +  a^d( r) r2  j  0_( r)  - 


tan 


-1  1 


L  P'  '  aod(r)  r  J 


+  V  1  +  a2(r)r2  [ep(  r)  -  tan-1  ^TTTf]  ~U~(7)  } 


(20) 


for  m  =  1  and  2. 


UQd( r)  =  local  speed  of  advance  in  off-design  condition 
( r)  =  local  speed  of  advance  at  design  condition 

The  values  of  ^ocj(  r)  ar|d  the  corresponding  aoc|(  r)  are  deter¬ 
mined  based  on  the  assumption  that  the  propeller  operates  in  all 
conditions,  design  and  off  design,  in  the  same  disturbed  fluid  field 
responding  with  different  propeller  performance  characteristics.  The 
change  of  value  of  advance  ratio  at  off-design  condition  is  assumed 
due  to  change  of  speed  of  advance  at  constant  RPM  according  to 

[UA(r>]od  [DA<r)]od  Jod 

[UA<r>]d  [0A<r)]d  Jd 
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from  which 

^7-tufl(r)]d 

where  U.(r)  and  U.( r)  are  the  local  speed  of  advance  and  its 

A  A 

volumetric  mean  value,  respectively.  Then  by  making  use  cf  the 
definition  of  the  advance  ratio 


a(r)  - 2— 

od  [uA(r)]c 


tUA(r)],  J 


the  advance  ratio  at  off-design  condition  is  given  by 

a(r)  =  a(r)  (22 

od  d  J0d 

The  normal  velocity  due  to  change  of  angle  of  attack  at  off-design 
condition  (AJ  =  -  J^)  can  now  be  determined.  Then  the  additional 

blade  loading  AL^°^  is  calculated  through  the  integral  equation  (5). 
od 

k.  Blade  Pressure  Distribution 

The  integral  equation  (5)  can  be  solved  successively  by  imposing 
boundary  conditions  given  in  Section  3  of  the  Analysis.  Then  the  span- 
wise  loading  distribution  is  determined  by 

r)  =  j  l/q)(  r,9  )  sin©  d0 
J  a  a  a 

o 

tt  n  m  a  x 

=  J  2  L(q'n)(r)  ©(n)  sin6adea  (2; 


For  the  unsteady  loading  due  to  wake  (q  f  o),  the  complete 
Birnbaum  series  is  used  for  the  chordwise  loading  distribution.  Thus 

rr  nmax 

L(q)(r)=^J  jL(q>1)(r)(^cos0a)  +  ^  X  L(  q>  n)(  r)s  i  n(  n- 1  )9& 

°  '  n=  2 

sin6o}deo=  L(q,,V)  +  >  L(q,2)(r)  (21*) 
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For  the  steady  state  loading  due  to  AJ  in  the  off-design  condition 
the  spanwise  distribution  is  obtained  by  the  first  term  of  Eq.  2k. 


In  the  case  of  NACA-a  mean  line  blade  sections,  with  "roof-top" 
blade  loading  distribution: 


L(0)(r)  - 


L(0,l)(r) 

/»  1+COS0 

L(0>'\r)  - 2! 

2(1-5) 


The  spanwise  loading  is  given  by 

L^(r)  =  L^°;  r)i  J  sinG^dB^  +  1 


0  s  0  <;  cos-,(  1  -  2a) 

a  ' 


cos-' ( 1-2a)  &  6  £  tt 

v  a 


TT 

f  ( 1+cosB  )s in0  dB  ^ 
"  '  cr  a  at: 


2(1-5)  X 


=  L(0jl)(r)  (5+0 


(25) 


where 

X  =  cos  '(1  -  2a) 

For  any  other  arbitrary  mean  line,  the  spanwise  loading  distri¬ 
bution  is  derived  from  Equation  (2k)  for  q  =  0 


tt  n  mO'x 

L^°,n)(  r)  =  1  J  -[\_( °>  1 )(  r)(  l+cos9^)  +  £  l/ °> n^(  r)s i n(  n- 1 ) 0^  sinG^jdB^ 

°  n=2 

=  L^'V)  +  L(0'2\r)  (26) 

Whenever  the  ^optangent  Birnbaum  mode  is  used,  the  Van  Dyke  or 
L i gh thill  correction  is  applied  to  the  chordwise  pressure  distribution 
to  remove  the  physically  unrealizable  leading  edge  singularity. 

Their  simple  rule  for  obtaining  a  uniformly  valid  solution  is 


C 

P 


S 

S  + 


C 


P 


(27) 
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where 

C  =  pressure  coefficient  according  to  first-order  thin 
P 

a i rfoi 1  theory 

S  =  chordwise  location  from  the  leading  edge  in  fraction  of 
chord  length 

pQ  =  leading  edge  radius  of  the  profile  in  fraction  of  chord 
length 


The  modified  chordwise  loading  distribution  is  then 


q\  r,  0  )  =  -  (■ 
'  *  or  tt  V 


1-COS0 


1+p  -cos 
o  a 


^-){L(^’)(r)cot^+  _^L(q'">(r)  sin(n-l)ej 

(28) 


n=2 


Since  the  pressure  near  the  leading  edge  is  governed  by  the 
first  term,  the  location  of  the  maximum  pressure  near  the  leading  edge 
is  determined  by 


d  (  Sin  ea 

d0  v  T  +  p  -  cose  J  ® 
a  o  a 

which  yields 

cos  0  =  — - - 

“  ,  +  qo 

This  is  the  choruwise  location  of  a  point  where  the  maximum  pressure 
near  the  leading  edge  will  occur  after  the  application  of  the  leading 
edge  correction,  provided  pQ  <  .015,  otherwise  the  sine  series  will 
influence  the  location  of  maximum  pressure  near  the  leading  edge. 
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ITERATIVE  PROCEDURE 

As  was  noted  earlier,  a  propeller  operating  behind  a  hull  is  subjected 
to  the  effect  of  a  hull  wake  of  varying  intensity  (weak  to  strong)  and  it  is 
also  influenced  by  its  self-induction.  The  former  can  be  taken  into  account 
by  a  harmonic  analysis  of  the  wake,  as  has  been  done  before;  the  latter,  self- 
induction,  however,  depends  on  a  priori  knowledge  of  the  propeller  loading 
which  requires  that  the  reference  surface  be  established  in  advance. 

The  harmonic  content  of  the  wake  based  on  the  local  speed  of  advance 
is  used  as  the  inflow  to  the  propeller.  The  wake  intensities  may  be  quite 
strong  and  even  of  the  same  order  of  magnitude  as  the  free  stream  velocity 
(ship  speed).  These  velocities  cannot  be  considered  as  perturbations  and, 
thus,  a  strictly  linearized  theory  is  no  longer  valid.  Second-order  terms 
cannot  be  ignored  as  small  quantities.  However,  second  order  theory  (nonlinear 
theory)  cannot  be  easily  developed  so  as  to  lead  to  a  manageable  numerical  ap¬ 
proach.  An  iterative  procedure  is  established  based  on  a  combination  of  the 
nonlinear  form  of  the  Bernoulli  equation  with  a  kinematic  boundary  condition 
existing  at  the  designed  operational  condition  for  the  selection  of  the  refer¬ 
ence  helicoidal  surface  around  which  a  perturbation  theory  will  be  developed. 

A  flow  field  closer  to  the  propeller  operating  condition  is  achieved  and  thus 
the  requirement  of  the  small  perturbation  quantities  will  be  reinforced.  This 
approach  is  applicable  to  moderately  to  heavily  loaded  propellers  immersed 
in  a  strong  wake. 


1)  The  Bernoulli  Equation 

For  a  moderately  loaded  propeller  operating  in  a  steady  inflow  field 
whose  velocity  components  in  the  axial,  tangential  and  radial  directions  are 
given  by 


d# 

U.  (  r)  +  - — - 
A  ox 


ril  +  Vy 


r  dcp 


W 


r 


respectively,  the  Bernoulli  equation  for  a  point  (x,r,«p)  in  the  field  becomes 

P0(*,r,cp)  i-  /  2  ,  di  2  d*  2- 

~f — +2i  (v  zr)  +(rW+V  7  WJ  AV  or)  ] 

o  , 

“  .  1  ..s  .  / 


57*2  L"«r*v  <29> 
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where 

p  is  the  mean  local  pressure 
ro  r 

=  local  speed  of  advance 

$  =  the  flow  velocity  potential  function 

and  VT  and  -  mean  wake  velocity  in  the  tangential  and  radial  directions, 
respectively. 

The  rest  of  the  symbols  are  as  defined  previously. 

It  is  assumed  that  for  mean  ship  wake  flow 

VT  =  0  (always  true  for  a  single  screw) 

a*0 

and  =  0  and  =  0  (no  radial  variation  is  taken  into  account) 

and  reference  pressure  p*  =  0. 


Thus,  Eq.(29)  yields 

64  . , ,64  .  2 


d§  64  N 2-,  -  .64  .,.64  2-,  p 

Lua  izr+^r)  j+L“rQ  7  w  +  ^7  %r)  ] = ~  77  (x’r^> 


which  includes  the  second  order  terms,  or 

r„  VV]d§o  "  \f\ &Vn  1  s$0  p0 ,  . 

LV  fc-yJ  sr  “L^  ’  i\7  arvj  7  37-  =  ”  77  (x- r>q,) 


(30) 


For  a  coordinate  system  fixed  in  space  and  for  a  right-handed  propeller  rotating 
with  angular  velocity  Q, 


cp  =  -Qt  . 

Thus,  the  above  equation  becomes 

,  ,64  .,64 


1  . 

I  1  o 


64 


lua  +  ixzrJl  or +  Lrn  ■  2\7  tq  _  77  (x>r>c;) 


Lett i ng 


.  64 

A  ( x ,  r ,  cp)  =UA  +  2  [oTJ 


d  4  _ 


B(x’r’^  rQ  Lr£i  -1(7^, 


dtp  /j 
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Eq .  (3  1)  i s  wr i t ten  as 

5$  d$  p 

A  dT  +  B  FT  "  ^  (x’r;t)  =  ^,(x*r;t) 

where 

PQ(X»  r»t) 

ij|(x,r;t)  =  -  -  =  acceleration  potential  function. 


(33) 


The  formal  solution  of  the  differential  Eq.(33)  is  obtained  by  the  method 
of  character i s t ics  provided  the  coefficients  A  and  B  are  known  constants. 

This  present  scheme  is  also  used  as  a  means  to  introduce  the  propeller 
induction.  Identifying  the  propeller-induced  velocity  in  the  axial  and 
tangential  directions  by 

dT=  ui  (x'rlt) 

and 

1  °  _  ,  ,x 
-7  W~  vi(x’r;t) 


respectively,  the  coefficients  of  differential  equation  (33)  become: 
A(x,r;t)  =  UA  +  |  u. 

B(x’r;t>  =  rD  Lr 0  +  T  vi] 


(34) 


These  induced  velocities  can  now  be  approximately  evaluated  by  taking  one- 
half  of  the  corresponding  values  at  infinity  (well-known  results  from  the 
momentum  theory).  The  steady  induced  velocities  in  the  far  field  can  easily 
be  calculated  once  the  propeller  loading  is  known.  In  fact,  in  Reference  7 
it  is  shown  that  the  steady  axial  velocity  at  x~* 00  in  uniform  flow  is  given  by 


-  N  ma?  ?A^(0)L(0’S)(r) 
"  ’  2npfUro  -t,  7 


(35) 


where  all  the  symbols  are  as  previously  defined  and  linear  dimensions  are  non- 
dimensional  i  zed  with  respect  to  propeller  radius. 


It  can  be  further  shown 


7 


that 


the  steady  tangential 


velocity  in  the  far 
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field  can  be  expressed  in  terms  of  the  axial  induced  velocity  by 

v. (x, r ; t)  =  -  —  u.(x,r;t) 
i  ar  i 

thus  the  coefficients  ( Eq . 3^+)  become: 


A(x, r;t)  =  ~  u.(x,r:t) 


,  *  1  .  „  i  U|(x,r;t)~l 

l(x,r;t)  =^Lrn  -  2  - - J 


(36) 


which  incorporate  the  propeller  speed  of  advance  and  the  propeller  induction, 
both  of  which  are  considered  to  vary  radially. 

As  will  be  seen  later,  once  the  potential  function  lQ(x,r;t)  is 
determined  at  a  stage  of  the  development,  the  blade  loadings  can  be  evaluated 
and  hence  the  propeller  induction  effects  as  well  as  any  other  flow  charac¬ 
teristics  can  be  calculated.  The  coefficients  given  by  Eq.(36)  are  considered 
to  be  known  at  that  stage,  and  hence  the  differential  equation 

9  ®  9  § 

A(x, r ;t)  +  B(x,r;t)  =  f(x,r;t)  (37) 

(see  Eq.33)  can  now  be  solved  by  the  method  of  characteristics  provided  that 
A  and  B  are  determined  at  x=0(at  the  propeller  plane) , i .e. , they  are  speci¬ 
fiable  constants.  The  characteristics  are  deduced  from 


dx  _  dt 
A  =  B 

which,  upon  integration,  yields 
Bx  -  At  =  constant 


Lett i ng 


5  =  Bx  -  At 
T1  =  x 

then 


(38) 


^  .  3  5  9  9 

37  "  S<  3T  +  \  5rT  =  8  31  +  StT 

3  s  9  y,  9  .9 

ot  "  5t  31  +  \  3ri  =  'A  3? 

Multiplying  the  above  relations  by  A  and  B,  respectively,  and  substituting 
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the  results  into  Eq.(37)  yields 
oi  . 

A  oTf  =  ^T'>r‘>  i  (B1K)3 

which  upon  integration  becomes 

*0(x,r;t)  =  ^  J  iCf.r;  £(Bt  -  £)  ]dT 
—  00 

Making  use  of  the  relations  given  by  Eq.(38),  the  above  solution  can  be 
written  as 

$Q(x,r;t)  =  ^  J  VCt , r ; t  -  |  (x-t) ]dT  (39) 

—  00 


This  is  of  the  same  structure  as  the  existing  relation  between  the  velocity 
and  acceleration  potentials  for  the  case  of  a  lightly  loaded  propeller,  which 
was  the  subject  of  the  previous  study. ^ 


The  above  relation  indicates  that  the  solution  $  (x,r;t)  for  the 
velocity  potential  is  obtained  through  the  so-called  "History  Integral" 
of  the  corresponding  acceleration  potential  where  the  time  is  shifted  by  the 

Q 

amount  —  (x-t)  .  Thus,  the  term  involving  the  "time  element"  becomes 


,  u.(x,r;t) 

rQrn-2  a  r 

rLuA4  Ui(x’r;t)] 


(x-t) 


m 


where  all  the  linear  quantities  are  expressed  in  a  non-dimensional  form  in 


terms  of  the  propeller  radius  r 


Thus,  the  "inverse  advance  ratio1 


r  .  .  v  2 

ar  =  —  which  was  introduced  in  the  study  of  the  lightly  loaded  propeller 

is  defined  in  the  present  study  by 


a  r  = 
e 


Ur  r-i^lL 
o  2  ar 


Vi  ui 


(41) 


which  will  be  called  an  "effective  inverse  advance  ratio’.1 


It  is  shown,  therefore,  that  the  present  problem  with  propeller 
exposed  to  the  advance  speed  and  influenced  by  its  self-induction  admits  of 
the  same  formal  solution  as  in  reference  2.  The  same  mechanics  can  be 
utilized  as  before,  provided  the  coefficients  of  the  differential  Eq.(33) 
are  known.  This  can  be  achieved  by  introducing  an  iterative  procedure. 
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2)  Iterations 

At  first  it  is  assumed  that  the  induction  effects  are  omitted  (i.e., 
u(*^=0  ,  vj^=0),  so  that  the  coefficients  A(x,r;t)  and  B(x,r;t)  become 

A(l,(0.r;t)  -  Uft(r) 

B(l>  -  I 


both  known  quantities,  in  which  case  the  solution  yields 


S(1)(x,r;t)  =  J  >  r ;  t  -  — y 


dT 


A  -® 


which  is  identical  to  that  of  a  "lightly  loaded"  propeller  except  that  the 
present  solution  incorporates  a  radial  variation  in  the  inflow  field.  The 
variation  does  not  introduce  any  additional  complication  since  the  problem 
is  solved  by  dividing  the  lifting  surface  into  radial  strips  with  specified 
inflow  field  at  each  radius. 

Knowing  the  velocity  potential  4^(x,r;t)  and  the  "effective  in¬ 
verse  advance  ratio"  a^  =  S-  =  (n  =  the  propeller  RPS,  and  a^  =  a), 


e  UA  UA  ‘  '  '  e 

through  the  ex i s t i ng 'machi ne ry 1 (ana lys i s  and  program)  the  steady  state  blade 

loading  L/?\,n^(r)  is  determined  (subscript  refers  to  the  order  of 
'  (2) 

iteration).  Then  the  induction  effects  u)  ;(0,r;t)  are  calculated  through 
Eq.(35).  A  new  set  of  coefficients  A^(0,r;t)  and  B^(0,r;t)  are  determined 
through  Eq.  (  36)  . 

(2) 

With  these  new  coefficients,  a  new  solution  4^  '(x,r;t)  is  established 


>(2) 


(2) 


(x,r;t)  -  t(2K’.r;t  -  J^yy  (x-t);b 


With  "effective  inverse  advance  ratio" 

u(2) 

(2i  “r°r  •  *  Tjnr 

the  new  loading  ^ ( r)  and  inductions  ul^  are  calculated,  and  thus 

(2)  ■ 
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establish  the  new  coefficients  A^(0,r;t)  and  8^(0,r;t).  At  that 
stage  the  new  iteration  begins. 

The  results  of  the  jth  iteration  are  presented  in  the  following 


table: 

Coefficients  A^  =  j  uj"*^ 


r  r 
o 


1 

2 


u 


(J) 


e 


J _ 

Qr  r 
o 


"Effect i ve 
i nverse 
advance 
coefficient" 


(j) 


...  Orr  -  y  a  ^'"^r 
,(j)=  _  o  2  e 


'•liv  +  i  “!j)’ 


<j) 


U  . 
I 


rl 


Solution  of  $  '  (x,r;t)  =  —  Q.y  J  ,  r ;  t  -  --Q-y  (x-T)^dT 

d . e .  A  A 


Steady- 
1 nduct ion 


(j) 


-  N 


2npfUr  - 
f  o  n=l 


max  n  (0)|_|  ■  ’  .)  (  r) 

t,  - — Uli) - 


This  iterative  procedure  will  continue  until  the  potential  function  £  (x,r;t) 

or,  more  conveniently,  the  "effective  inverse  advance  coefficient"  ,  does 

not  vary  in  two  consecutive  iterations.  With  the  final  value  of  a^-^ ,  the 
reference  helicoidal  surface  is  established  and  this  will  be  used  as  the  basic 
information  for  the  calculations  of  the  steady  and  unsteady  blade  loading  and 
hydrodynamic  forces  and  moments. 


It  is  clear  that  the  number  of  required  iterations  depend  to  a 

great  extent  on  how  close  the  initial  value  of  the  "effective  inverse  advance 

coefficient"  a^  is  to  its  final  value.  A  simplified  procedure  is  given  below, 
e 

It  is  known  that  the  pitch  of  a  well  designed  1  Wake-adapted  propeller" 
nearly  satisfies  the  following  "kinematic  condition" 
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+  u.  =  (Or  +  v.)tan0p 


=  (Ur  -  — )  tanG 
'  ar  p 


from  which  the  axial  induction  is 

u.  =  (Qr  tanG  -  U.)  - ^-r — -r-  (42) 

i  '  p  A'  a  r  +  tanG  v  ' 

P 

where  9^  is  the  propeller  pitch  angle. 

From  Eq .(41)  the  axial  induction  is  determined  as 

u  =  2(rQ  -  arUn)  - - r  (43) 

1  1  +  a  r 

which  is  considered  to  be  the  result  of  the  "dynamic  condition"  since 
Eq.(4l)  is  derived  from  the  Bernoulli  equation. 

Combining  the  "kinematic"  (Eq.42)  and  "dynamic"'  (Eq.43)  conditions, 

yields 

E  a2r2-  2D  ar  ••  E  =  0  (44) 


whose  solution  leads  to 


(45) 


where 

D  =  Ur  -  U.  tan9 
A  p 

E  =  Ur  tan9  +  U. 

p  A 

With  the  initial  value  of  a  obtained  from  Eq .  (45) ,  the  iterative  procedure, 
described  above,  starts.  It  is  seen  from  the  numerical  calculations  that 
this  approach  has  reduced  the  number  of  iterations  cons i derab 1 y, by  50  percent. 

Once  the  reference  helicoidal  surface  is  established  and  the  "effective 
inverse  advance  coefficient"  ag  is  determined,  the  numerical  procedure 
becomes  identical  to  that  developed  previously  in  Reference  2  for  the  steady 
and  unsteady  flow  conditions. 
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PRESSURES,  HYDRODYNAMIC  FORCES  AND  MOMENTS, 

AND  BLADE  BENDING  MOMENTS 

1)  Blade  Pressure  Distribution  on  Each  Blade  Face 

In  the  preceding  sections  the  pressure  difference  Ap^(r),  at  a  given 
radial  position  r  and  frequency  q  ,  between  back  (suction)  and  front  (pres 
sure)  faces  of  the  blade  surface  is  determined.  On  the  suction  side  of  the 
pressure  due  to  loading  is  +Ap^  ( r) /2 ;  on  the  pressure  side  it  is  -Ap^q\r)/2 

In  addition,  a  non-lifting  pressure  P^  is  generated  due  to  the  sym¬ 
metrical  "flow  distortion"  thickness  effect  which  will  be  present  only  in  the 
steady  state  since  the  blade  is  considered  to  be  rigid.  PT  is  derived  in 
Appendix  C  by  means  of  the  "thin  body"  approximation  (see  Reference  8), 


The  instantaneous  pressure  is  the  sum  of  the  blade  pressures  due  to  all 
frequencies  contributing  significantly.  On  the  pressure  or  suction  face, 
when  the  blade  swings  around  its  shaft  in  the  clockwise  direction  from  its  up¬ 
right  position  (12M),  it  is 


P^s(r> 


r!qQ  =  Z  lp(q)s(r)lcos(qe-cp  ) 
q=o  P’  q 


where 


cp  =  phase  angle  =  tan 


(ptqK 

_pJsj_lm 

(PK) 

'  p.s'Re 


(the  subscripts  Re  and  Im  indicate  the  real  and  imaginary  parts)  and  0 
is  blade  angular  position,  positive  in  the  counterclockwise  direction 


as 


is  cpq. 


2)  Propeller-Generated  Forces  and  Moments 

The  principal  components  of  the  hydrodynamic  forces  and  moments  are 
shown  in  Figure  2  with  the  sign  convention  adopted.  The  total  forces  at  fre¬ 
quency  £N(£=0, 1,2,. . .)  induced  by  an  N-bladed  propeller  are  determined  as 
(see  Reference  2)  : 
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F  ■ 


F  = 


Re  {NiQ  e,SlNnt  }  L{*N)(r)cos6  (r)dr) 
o  P 


H>  .1  -  _  . 

Ke  -e,tKnt  /  [  fL(«-l.n)(r)A{n)(.e»-)  +LU«+l.n)(  )A(S)(6r)3 


S!n6p(r)dr> 


Fz  =  Re  {^-ei£Nnt  }  $  iL^'1  •">  (r)A("}  (-9^) -  L(lN+1 '"}  (r)A("}  {0^)1 


o 

n=  I 


sin0  (r)dr) 
P 


(^a) 


The  moments  are  determined  by: 

=  -Re  {N  r^  el£Nat  J  L^N^(r)sin9  (r)rdrl 


Q  =Re{^°  e!^t  /  {  J  [LUN-1,n)(r)A(n)(_0r)+L(^l,n)(r)A(n)(9r)]cos,rr) 
°  n=l 

+  l  [L(£N"1’n)(r)A1(n)(-0')-L(lN+,’n)(r)41(n)(e^)]Cie')Sin0  (r)tane  (r)}rdr} 

lb  b  b  p  p 


n=  1 
-Nr 


/  {  l  [L(“",'n)(r)Aln)(-e^)-L(£“+l’n)(r)A(n)(8^)]co50t,(r) 


Q  *  , 
n- 1 


l  [LUN'1’n)(r)A1(n)(-S^H(tN+,’n)(r)A1(n)(ep](i9^)  surf  (r)tan0  (r)}rdr} 


n=l 


P  P 

(47b) 


where  (z)  and  (z)  are  as  defined  in  Appendix  A. 

It  is  seen  from  Eqs.  (47a, b)  that  t fie  transverse  bearing  forces 
and  bending  moments  are  evaluated  from  propeller  loading  components  L^’^Jr) 
associated  with  wake  harnionics  at  frequencies  adjacent  to  blade  frequency, 
i.e.  at  q  -  IN  -  1 ,  whereas  the  thrust  and  torque  are  determined  from  the 
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loading  L^(r)  at  blade  frequency  q  =  AN  .  At  X=0  (steady-state),  the 
mean  transverse  forces  and  tending  moments  are  determined  at  shaft  fre¬ 
quency.  Thus: 


Re  {^/I  l  LO’n)(r)A(n)(6[))sin0p(r)dr} 

O  -  ! 

n=l 

Re  {^f/1  l  L(,'n)(r)A(n)(6[)]sin9p(r)dr} 
o 

n~  • 

6  =  Re  {-— •  f  i[  l  L(l>n)(r)A(n)(erb)  cos0p(r)  - 
°  n=l. 

[  l  L(1  ,n^  (r)A^  (6^)3  (i6b)sin0p(r)tan0p(r)}rdr} 
n  =  1 

Qz»  Re  }  U  t  L(l’")(r)AI")(e')  cosMO  - 

o  ~  . 

n-I 

[  l  L^1  ,n^  (r)A^  (6^)]  (i6b)sin0p(r)tan0p(r)}rdr}  (*+ 7c ) 

n-l 


All  forces  and  moments  can  be  written  in  the  form 
Re[C^e  9q  e',q°]  =  C^cos(-qO  +  cpq)  =  C  tq^cos  (qO  -  T>q) 


where  0  is  blade  angular  position,  red  live  Ik  the  cov.ntcr-cloctavise 
direction  from  zero  at  the  upright  position  (12M),  q  is  order  of  shaft 
frequency,  C^q^  is  magnitude  of  force  or  moment,  and  Cfq  is  phase  angle 
(electrical)  determined  as  the  angle  whose  tangent  is  the  imaginary  part 
(sine  part)  over  the  real  part  (cosine  part).  (The  program  output  gives 
C  ^  and  cpq.  ) 
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When  9^  i s  pos i t i ve, 
the  upright  position,  i.e. 


the  peak  (or  trough)  of 
leads:  when  cp  is  negative, 


is  to  the  left  of 
the  peak  (or  trough) 


is  to  the  right  of  the  upright  position  and  lags  that  of  input. 


In  Reference  2,  as  a  first  attempt  at  estimating  the  frictional  con¬ 
tribution  to  thr  ust  and  torque  Q.  ,  use  was  made  of  the  Prandt  1 -Schlicht  i ng 

formula  for  the  friction  coefficient  for  one  side  of  a  smooth  flat  plate  as 

9 

was  suggested  by  Hoerner.  'I  he  present  study  follows  Reference  2  in  this 
matter  with  this  exception:  in  the  formulas  for  friction  coefficient  and 
frictional  thrust  and  torque,  the  inverse  advance  ratio  is  no  longer 
a  =  Qro/U,  based  on  forward  speed  of  the  ship  and  constant  over  the  propel¬ 
ler  span,  but  a^(r)  =  Or/U^(r),  based  on  local  speed  of  advance  and  varying 
with  radial  position. 


3)  Blade  Bending  Moments 


The  blade  bending  moment  about  the  face  pitch  line  at  any  radius  r. 
of  a  blade  is  calculated  from  the  chordwise  integrated  loading  (spanwise 
component)  L^(r)  at  any  shaft  frequency  q  as 


W«nt  =  r3  eiqnt 

D  O 


(48) 


The  positive  blade  bending  moment  about  the  face  pitch  line  is  that  which 
puts  the  face  of  the  blade  in  compression. 

The  instantaneous  blade  bending  moment  distribution  when  the  propeller 
swings  around  its  shaft  in  the  clockwise  direction  is 


where  cp^  is  the  phase  angle  (electrical). 

It  should  be  noted  that  in  the  program  the  value  of  r.  of  Eq.(48)  is 
limited  to  any  of  the  midpoints  of  the  radial  strips  into  which  the  blade 
span  is  divided,  at  which  points  the  pitch  angles,  as  well  as  other  geomet¬ 
rical  characteristics,  are  given  as  input.  The  bending  moment  at  any  other 
radial  position  can  be  obtained  by  interpolation  or  extrapolation. 
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NUMERICAL  RESULTS 

The  theoretical  approach  developed  in  the  preceding  sections  for  the 
evaluation  of  the  blade  pressure  distributions  due  to  loading  and  thickness 
effects  and  of  the  resulting  forces  and  moments  has  been  adapted  to  a  high¬ 
speed  digital  computer  (CDC-6600  or  Cyber  176). 

The  numerical  procedure  has  been  applied  to  propellers  of  different 
loading  from  light  to  moderate  to  heavy: 

a)  the  DTNSRDC  Propeller  4 1 1 8  tested  at  DTNSRDc'0’^  in  screen  wakes 
(3-cycle  and  4-cycle  screens) 

b)  the  Sharp  Propeller  V-3275  in  the  wake  of  ship  model  4986  (S.S. 
MICHIGAN),  the  wake  survey  and  propeller  geometry  supplied  by  the  U.S.  Mari¬ 
time  Administration 

c)  the  NSMB  Propeller  4930  in  the  wake  of  ship  model  4705,  the  wake 
survey  and  propeller  geometry  supplied  by  Lloyd's  Register  of  Shipping 

The  particulars  are  listed  in  Table  1. 


TABLE  1 

Propeller  Designation 

DTNSRDC  4118 

SHARP  V-3275 

NSMB  4930 

Expanded  Area  Ratio,  EAR 

0.60 

0.564 

0.810 

Pitch-Diameter  Ratio, 

P/D  at  0.7  Radius 

1.077 

1.080 

0.757 

Skew,  deg 

0 

11.5 

30.0 

Diameter,  D,  ft 

1 .0 

22.5 

6.837 

No.  of  B  lades,  N 

3 

5 

5 

NACA  a  Meanline  Section 

a  =  0.8 

a  =  0.4 

a  =  0.7 

RPM  (=  60n) 

900 

106 

213.6 

Advance  Ratio,  J  =  U/nD 

0.831 

0.993 

0.746 

Free  Stream  Velocity,  U, 
ft/sec 

12.465 

39.50 

18.158 

Speed  of  Advance,  U.,  at 
0.7  Radius , ft/sec 

12.81 

33.18 

9.02 
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Relevant  geometric  characterist ics  of  the  3  propellers,  namely,  the  ratios 
of  maximum  camber  to  chord  length,  m^/c,  of  maximum  thickness  to  chord, 
t  /c,  of  leading  edge  radius  to  chord,  pQ,  of  chord/diameter,  c/D,  and  of 
pitch  to  diameter,  P/D,  are  given  in  Tables  2-4. 


TABLE  2. 

DTN5RDC 

PROPELLER  4118 

Radius 

m  /c 

X 

t  /c 
0 

Po 

c/D 

P/D 

0.25 

0.0228 

0.090 

0.00525 

0.347 

1.085 

0.35 

0.0231 

0.068 

0.00290 

0.386 

1.083 

0.45 

0.0224 

0.052 

0.00170 

0.425 

1.081 

0.55 

0.0212 

0.040 

0.00100 

0.454 

1.079 

0.65 

0.0203 

0.031 

0.00060 

0.464 

1.077 

0.75 

0.0198 

0.024 

0.00035 

0.452 

1.076 

0.85 

0.0189 

0.018 

0.00025 

0.405 

1.074 

0.95 

0.0174 

0.016 

0.00020 

0.278 

1.072 

TABLE  3: 

SHARP  PROPELLER  V- 3275 

0.25 

0.0553 

0.232 

0.059 

0.166 

1 .000 

0.35 

0.0445 

0.166 

0.030 

0.200 

1.032 

0.45 

0.0374 

0.122 

0.016 

0.231 

1.055 

0.55 

0.0321 

0.094 

0.0097 

0.251 

I.083 

0.65 

0.0268 

0.073 

0.0059 

0.260 

1.080 

0.75 

0.0218 

0.056 

0.0035 

0.257 

1.079 

0.85 

0.0174 

0.042 

0.0019 

0.236 

1.073 

0.95 

0.0138 

0.031 

0.0011 

0.168 

1.057 

TABLE  4. 

NSMB  PROPELLER  4930 

0.232 

0.0304 

0.168 

0.04087 

0.234 

0.729 

0.335 

0.0288 

0.124 

0.02238 

0.268 

0.750 

0.437 

0.0279 

0.088 

0.01210 

0.308 

0.754 

0.539 

0.0258 

0.065 

0.00680 

0.343 

0.758 

0.642 

0.0217 

0.049 

0.00309 

0.364 

0.759 

0.744 

0.0178 

0.039 

0.00183 

0.366 

0.754 

0.847 

0.0138 

0.030 

0.00162 

0.343 

0.741 

0.949 

0.0101 

0.023 

0.00175 

0.264 

0.724 
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Mean  wake  information  is  presented  in  Table  5  as  the  ratio  of  speed 
of  advance,  1)^,  to  the  free  stream  velocity,  U,  at  8  equidistant  radial 
positions  between  propeller  hub  and  tip  (see  Tables  2-4). 


TABLE  5.  UA/U  -  RATIO  OF  SPEED  OF  ADVANCE  TO  FREE  STREAM  VELOCITY 


Propeller  DTNSRDC  4118 

SHARP  V-3275 

NSMB  4930 

Wake 

Screen 

Mode  1  4986 

Model  4705 

U,  ft/sec 

12.465 

39.50 

18.158 

Radial  Position 

Vu 

UA/U 

vu 

1 

0.970 

0.599 

0.207 

2 

0.969 

0.678 

0.285 

3 

0.984 

0.747 

0.35A 

4 

1.009 

0.800 

0.407 

5 

1.026 

0.829 

0.452 

6 

1.028 

0.844 

0.527 

7 

1.026 

0.848 

0.655 

8 

1.025 

0.852 

0.791 

DTNSRDC  Propeller  41 18 

Results  of  calculat 

:ions  of  the 

hydrodynamic  forces  and 

moments  by 

the  present  approach,  as 

implemented 

by  the  PLEXVAN  program, 

are  presented 

in  Table  6  for  the  3-bladed  propeller  4118.  The  steady  and  blade-frequency 
K  coefficients  are  compared  there  with  available  experimental  values,'^  and 
with  values  obtained  through  the  approach  of  Reference  2  by  the  PPEXACT 
modified  code.' 


Since  the  publication  of  Reference  2  the  numerical  procedure  of  the 
PPEXACT  program  has  been  modified  for  greater  accuracy  of  the  chordwise 
d i s  t  r i bu t i ons . 
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TABLE  6 

DTNSRDC  PROPELLER  41 18 
IN  3-CYCLE  AND  4-CYCLE  SCREEN  WAKES 


CALCULATIONS 

EXPERIMENTAL 

PLEXVAN 

Program 

PPEXACT 

Program 

Measurements 

Coefficient 

Mag. 

Phase 

deg 

Mag. 

Phase 

deg 

Mag. 

Phasi 

deg 

Mean 

Rt 

0.166 

180 

0.166 

180 

0.154 

180 

0.0324 

0 

0.0323 

0 

0.0290 

0 

Blade 

Frequency 

RT 

0.0735 

-  22.9 

0.0742 

-  23.4 

0.0685 

-32 

O’ 

0.0126 

157-2 

0.0127 

156.6 

0.0102 

135 

Kr 

0.0175 

-  8.3 

0.0162 

-  8.9 

0.0139 

30 

y 

V 

0.0124 

-  24.5 

0.0114 

-  25.0 

0.0109 

30 

y 

KFz 

0.0175 

-  98.3 

0.0162 

-  98.9 

0-0131 

-67 

\ 

0.0124 

-114.5 

0.0114 

-115.0 

0.0114 

-67 

Comparisons  between  the  results  of  the  two  theoretical  approaches  for 
the  chordwise  and  spanwise  distributions  of  pressure  are  graphically  ex¬ 
hibited  in  Figures  2  to  4. 

It  is  seen  that  for  this  lightly  loaded  propeller  with  =  8k^./ttJ'^ 
approximately  0.6,  there  are  only  slight  differences  between  the  present 
approach  and  that  of  Reference  2,  as  was  expected. 
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Sharp  Propeller  V-3275 

Similar  calculations  have  been  performed  for  the  Sharp  5-bladed  pro¬ 
peller  V-3275.  The  results  are  shown  in  Table  7  and  graphically  in  Figures 
5-8.  Also  tabulated  are  some  results  obtained  from  Troost  B-5  propeller 
charts . 

TABLE  7.  SHARP  PROPELLER  V-3275  IN  THE  WAKE  OF  THE  S.S.  MICHIGAN 

CALCULATIONS 


PLEXVAN  Program 

PPEXACT 

Program 

TROOST  B-5  Prop 

Coeff i c i ent 

Mean 

Mag.  Phase 

deg 

Mag. 

Phase 

deg 

Mag. 

*T 

0.222 

1 80 

0.226 

180 

0.185 

O’ 

1^ 

0.0400 

0 

0.0408 

0 

0.0325 

i<F 

y 

0.0159 

180 

0.0171 

180 

\ 

0.0111 

180 

0.0118 

180 

Rpz 

0.0144 

0 

0.0159 

0 

^z 

Blade 

Frequency 

0.00983 

0 

0.0103 

0 

*T 

0.00635 

177.8 

O.OO638 

166.9 

~Ko 

0.00108 

-  2.1 

0.00106 

-  12.9 

kf 

y 

0.01265 

-140.8 

0.01171 

-143.1 

0.00842 

-153.1 

0.00777 

-152.5 

~Kpz 

0.00442 

-  75.2 

0.00437 

-  76.1 

~^z 

0.00274 

-  74.7 

0.00272 

-  73.2 

•A. 

“at  =  U^/nD  at  0.7  radius 


In  this  case,  C^.  =  8  Kj/nJ^2  is  approximately  0.8.  It  is  seen  that 
the  hydrodynamic  force  and  moment  results  of  the  present  analytical  approach 
differ  slightly  from  those  of  the  PPEXACT  modified  program.  The  blade  pres¬ 
sure  distribution  of  Figures  7  and  8  show  larger  differences  between  the 
two  methods. 
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NSMB  Propeller  4930 

A  comparison  of  results  of  the  PPEXACT  and  PLEXVAN  programs  for  this 
5-bladed  propeller  is  presented  in  Table  8  and  Figures  9  to  II.  The  compu¬ 
tations  were  performed  for  the  steady-state  case,  q=0,  only,  since  just  the 
mean  wake  velocity  was  supplied. 

Both  Davidson  Laboratory  sets  of  calculations  were  compared  with  the 
results  of  the  NSMB  lifting-surface  program  which  were  provided  in  the 
form  of  graphs  of  spanwise  load  distribution  and  chordwise  pressure  distri¬ 
bution.  The  values  of  mean  thrust,  K^,  and  torque,  K^,  coefficients  shown 
in  Table  8  for  the  NSMB  program  are  derived  from  the  spanwise  loads. 

This  propeller  with  a  Cy«4  is  more  heavily  loaded  than  the  other  two 
propellers  and  the  results  show  larger  differences  when  the  two  theoretical 
approaches  are  compared. 


TABLE  8.  NSMB  PROPELLER  4930  IN  THE  WAKE  OF  MODEL  4705 


Mean  Coefficient 

PLEXVAN 

CALCULATIONS 

PPEXACT 

NSMB 

TR00ST 
B-5  PR0I 

Rt 

-0.202 

-0.216 

-0.227 

-0.20 

0.0272 

0.0294 

0.0305 

CNJ 

o 

o 

at  J, 


U^/nD  at  0.7  radius. 
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CONCLUSIONS 

A  theory  and  corresponding  computer  program  have  been  developed  for 
a  marine  propeller  operating  in  non-uniform  inflow  by  taking  into  account 
the  radially  varying  mean  wake  and  mean  propeller  induction. 

The  helicoidal  reference  surface  of  radially  varying  pitch  is 
established  through  an  iterative  procedure  by  combining  the  non-linear 
form  of  the  Bernoulli  equation  with  an  appropriate  kinematic  condition  ex 
isting  at  the  design  stage. 

The  numerical  procedure  furnishes  information  on  the  blade  pressure 
distribution,  as  well  as  the  resulting  hydrodynamic  forces  and  moments  in 
the  steady  and  unsteady  flow  conditions.  The  approach  is  valid  for  mod¬ 
erately  to  heavily  loaded  propellers. 

It  is  seen  from  this  small  number  of  calculations  for  propellers 
with  thrust  coefficient  C^  varying  between  0.6  and  4.0  that  the  dif¬ 
ferences  between  steady  and  unsteady  forces  and  moments  calculated  by  the 
present  theory  (coded  as  PLEXVAN)  and  the  corresponding  values  determined 
by  the  previous  approach  (coded  as  PPEXACT)  are  larger  as  C-j.  increases. 
However,  the  differences  between  the  results  of  both  procedures  are  minor 
even  for  the  propeller  operating  under  more  heavily  loaded  conditions. 

The  PPEXACT  method  assumes  the  reference  surface  (along  which  the 
shed  vorticity  is  considered  to  be  convected)  to  be  a  pure  helicoid 
whose  local  pitch  is  fixed  by  the  joint  action  of  the  forward  ship  speed 
and  the  blade  tangential  velocity  at  any  radius.  The  difference  between 
the  ship  speed  and  the  radially  varying  wake  speed  is  taken  as  a  pertur¬ 
bation  from  this  surface.  On  the  other  hand,  in  the  present  (PLEXVAN) 
method  to  begin  with  the  reference  surface  (based  on  the  non-linear 
form  of  the  Bernoulli  equation)  takes  into  account  the  radially  varying 
mean  wake  and  mean  propeller  induction. 

It  would  appear  that  the  two  methods  yield  comparable  results  pro¬ 
vided  that  the  prescribed  procedures  are  closely  adhered  to  (e.g.,  in  the 
PPEXACT  method  the  hydrodynamic  pitch  tan  (U/rfl)  must  use  U  =ship  speed) 
The  advantage  of  the  PPEXACT  program  lies  in  its  considerable  saving  of 
computer  time. 
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If,  however,  greater  accuracy  in  the  evaluation  of  hydrodynamic 
forces  and  moments  is  required  and  more  precise  blade  pressure  is  needed, 
then  the  program  PLEXVAN  must  be  exercised. 

In  the  analysis  performed  in  Reference  12  to  evaluate  the  sensitivity 
of  various  wake  parameters,  the  hydrodynamic  pitch  angle  @  was  found 
to  be  the  most  critical  one.  A  set  of  calculations  for  propellers  with 
36°  skewness  indicates  that  a  decrease  in  @  of  10°  results  in  a  decre¬ 
ment  of  approximately  1 8%  in  K^.  (see  Figure  11  of  Reference  12).  On  the 
other  hand,  the  present  calculations  for  the  NSMB  4930  propeller  with  30° 
skew  indicate  that  with  a  decrease  in  3  of  4.5°,  there  is  a  decrement 
of  approximately  7%  in  K^.  The  trend  is  shown  in  both  cases  to  be  about 
the  same.  For  a  final  assessment,  more  systematic  calculations  must  be 
performed  and  compared  with  corresponding  measurements. 
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FIGURE  I.  DEFINITIONS  OF  ANGULAR  MEASURES 

NOTE:  THE  BLADE  REFERENCE  LINE  IS  THAT  CONNECT  ING  THE  SHAFT  CENTER  WITH 
THE  MIDPOINT  OF  THE  CHORD  AT  THE  HUB 
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STEADY-STATE  CHORDWISE  PRESSURE  DISTRIBUTION 
AT  0.65  RADIUS  ON  DTNSRDC  PROPELLER  4118  AT  DESIGN  J 
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FIG.3.  BLADE  FREQUENCY  SPANWISE  LOAD  DISTRIBUTION 
FOR  DTNSRDC  PROPELLER  4118 
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BLADE  FREQUENCY  CHORDWISE  PRESSURE  DISTRIBUTION  ON  SUCTION  AND  PRESSURE  SIDES 
FOR  PROPELLER  V3275  AT  r/ro=0.65 
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FIG.!  I.  STEADY  STATE  CHORD  WISE  BLADE  PRESSURE  DISTRIBUTION 
ON  THE  PRESSURE  AND  SUCTION  SIDES  AT  r/r0  =0.744 
FOR  PROP.  4930 
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APPENDIX  A 


Evaluation  of  the  cp^-  and  9  -  Integrals  of  the  Integral  Equation  (6) 

dcp 


tt 


1)  |W|y)  .1  J  }®r,vcos'' 


l(l)(»)  -if  (I -cost?)  e 'ycoscP 


=  JQ(y)  -  'J,  (y) 


I  ^  (y)  =  ~  J  ( 1  +2coscp)  e  1  yc0SCP  dcp  _  j  (y)  +  j 2 J  (y) 


I  (l">2)  (y)  -  i  J  costi-Dq.  e  'VC0S'P  dq>  ,  i”'1  j_  (y) 
o  m-1 

where  J^fy)  is  the  Bessel  function  of  the  first  kind  of  order  n  and 
argument  y 

(fi) 


rr  -izcosG 

2)  A^n;  (z)  =  J  S(n)  e  sinp  d6 

o 

a)  Birnbaut)  distribution 

(l  )  ,  ^  _  i  r™  __  9  _  -  i  zcosB 


au;(2) -U 


cot  -  e 


s  in0  d0  =  Jq(z)  -  i  Jf  (z) 


A^n>'  ^  (z)  =  —  f  s  in  (n-1 )  9  sin9e 

TT  J 


- 1 zcos  9 


d9 


nr* 


[j.  (z)  *  J.(z)] 


n-2  n 

b)  "Roof-top"  distribution  (a  mean  lines) 

(0(z)  =  '■C°S’,(l'2a") 
o 


A'  ' (z)  -  J 

o 

T! 


-  i  zcos  9  .  Q 

e  s  i  n9  d9 


cos  (1 -2a) 


(l+cosD)  -i zcos 9  .  .  .. 

2TiTi)  e  s,n6de 


„  -  i  z  l  _i.  .  1  f  i  2az 

'  )  2  2 ( 1 -a) zz  ,e  "  6 


j-ei2az  .  ei2zJ  | 


A(n>,)(z)  =  0 


(for  a® l,  A^^(z)  =  —  - n  *  ,  and  for  z=0,  A^(0)  =  l+a.) 
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Functions  Required  for  Evaluating  the  lntegrandof  the  Kernel  Function 
at  the  Singularity  (see  Reference  2)  and  the  Propel  ler-generated 
Moments  (see  text) 

1)  ~  f  e'^co5^  costp  dtp 

o 

>,0)(y)  -  -j[y»)  -  J2(y)]  +  u,<y) 
if2)(y)  =  [j„(y)  -  J2(y)]  +  iJ,(y) 
i,fe2,(y)  -V[-  d.(y)  (y)] 

m  m~2 

2)  Aj  (z)  =  J  S(n)  sinG  cosG  e  '2COS0  jq 

o 

a)  BirnbaLm  distribution 


A,0)(z)  =l[Jo(z)  '  J2(z)]  '  ?Vz) 


b)  "Roof-top  distribution  (a  mean  lines) 


A(i) _  -iz  n  ,  _i_ _ l_ 


I 


f  1  ~2a  _  j.)  i2az 
2  z  e 


+  (hi)*i2z]} 


(For  a  =  1,  Aj^(z)  =  ( 


cos  z  -  ii-n— i)  .) 


It  is  to  be  noted  that  the  values  for  negative  argument,  i*e*, 

I  ^  ^  (~y)  .  1  (-y)  .  A^(-z)  and  A^  (-z)  ,  are  the  conjugates  of  the 

values  given  above* 
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text 


APPENDIX  B 

Normal  Velocity  Due  to  “Nonplanar11  Blade  Thickness 
The  following  substitutions  are  made  in  Equations  (15)  and  (16)  of  the 


d0  =  G^sinG  d9 
o  b  a  a 


3y§,p>0o)  aA(P)  afT(p,ea) 
5f 


Gl s i nG 
b  a 


a 


where 


,  ,  »  im(9  -cp  +0  )  ®  ,  e\ 

1  =  1  S  e  °  °  n  f  .  .  iu(x-g) 

R  TT  m-_00  J  UK)m  6 


/  lm(|u|p)Km(|u|r)  for  p<r 

(IK)  =  { 

m  ll(lulr) 


du 


a 

On 1 


,(|u|  r) Km ( | u |  p)  for  r  <  p 

T===-  (v«  5«  -  7 

^1+a*(r)  r2  '  °/ 


..  V  im^n  fN  for  m=£N,  J-0,±1,±2,  ... 

Since  e  n  =  l  ,  Equation  (16) 

n~l  *  0  otherwi se 


becomes 


(°) 

VT  (0 

uTTT 


iNr 


2^3UA(r)Vi+a2(r)  r2  ^  °p 


^  rr„  - -  Sf  i m ( 0  -cp  ) 

5  jl  UA(p)Vl+a2(p)p3  _!e  o  o' 


j\aA( r) u  +  ^rjo^m®1  (x'§)ududpd0a  (B-l) 


With  the  trigonometric  transformations  for  x,  cpQ,  and  9q  given  in 
the  text,  the  series  is  brought  to 


B  1 


R-2063 


JL  I  J  UA(c)Vl+a=(p)p^ 

m=i'N  °  P 


dfT 

Se7 


■im(or-ap) 


+ 


iu(ar/aA(r)-ap/aA(p)) 


'  (u-aA(p)m)ePcosea/aA(p) 


■i,"»(r',)'^Va»|rl(,dpd»a 


( B  — 2) 


Let 


TT  df 

A((u-a  (p)m)eP/a  (p))=  J 

o  a 


'  ((u-aA(p)m)6p/aA(p))cos0a 


and  applying  the  generalized  lift  operator  (see  Appendix  A) 


V 


w 


_ i_Nr _ 

2n2uA(r)^l+a^(  r)  r2 


J  UA(p)^l+aA(p)p3 


•£  (r)u+a)(IK)  iu(0r/aA(r,  -ap/aft  (p) ) 

X=-oo  -oo  r 


'XN 


a(  (u-aA(p)  XN)  0P/aA  ( p)  j  1^  (( aA  ( r)  Xn-u)  d^/a^  (  r))dudp  (B  -3) 


It  can  readily  be  seen  that  folding  the  doubly  infinite  X-series  and 
u-integration  yields  Equation  (17)  of  the  text. 
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APPENDIX  C 


Thickness  Effects  (Non-lifting) 

The  velocity  potential  due  to  source-S-ink  distributions  approximating 
the  symmetrical  thickness  distribution  of  an  N-bladed  propeller  is  given  by 


*(x,  r,cpQ;0) 


1 

An 


N 

£ 

n=1 


M(C,s) 

R 


dS 


^  f  H(p,s)  ^l+a2(p) p2 
JJ  R  a(p)p 


p6bsin9adGadp 


(C-l) 


where  the  source  strength  is 


, - - - r  SfT(P.s)  of  (§,p,0  ) 

M(p,s)  =  2UA(p)^l+a‘'(p)p2  - -  =  2UA(p)  ■ 


the  Descartes  distance  R  is 

2  ^  /  -j£ 

.  (x-§)  +  r“  +  p2  -  2rp  cos\6  -cp  +§  )  | 

l\  /  r  \  o  To  n/_l 

s  is  measured  along  the  chord  of  the  blade  section  and  the  other  symbols 
are  as  defined  in  the  text. 

From  Bernoulli's  equation,  the  linearized  pressure  P^  is 


Mx’r’Cpo>  =  -pfUA(r)  +  a(r>  ^ 


Hence 


j^.  /.xdfT(5’p,e0) ;yd 

op 


PT(x’r’Cpo)  =  — IJUA  ^  p)  ~  a  rpT^C"  LW  +  3  ^  r>  j  0b  ^ 1  +a'  < p >  p"  sinSodVp 

(C-2) 


If  use  is  made  of  the  expansion 


1  1  .*  im(90-<p  +9  )  • 

IEi  e  n'  J(1K) 

m=-co  -co 


ei(x-5)kdk 


where 


i  I  (  I  k I  p)  K  ( I  k I  r)  for  p  <  r 

(lK)m  =  m 

I  1 m(  t  r) Km( 1 k I p)  for  r  <  p 


Cl 


t 
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a  A  v  a  1  _i  ?  _  im(e0-'?  +?0) 


&+a<r> 


ocp  J  R  n 
*  o  m=  -  ; 


J(k-a(  r)m)  (  IK) 


(x-s)  k 


S  i  nee 


(C-3) 


%  =  0Q/a(p)  =  (aP  -  0P  cos0a)/a(p) 


dfT(5*P-60) 


.  .  df  (p,e  ) 

a ( p)  txh’  a7 


6b  sin9o 


(C-4) 


N  ;m0  (N  for  m  =  IN  ,  £  =  0,  =1,  =2,  ... 

X)  e  n  =  (  „  „  ,  „ 

n=  1  0  for  m  /  a 


(C-5) 


On  substituting  Eqs.(C-3)  -  (C-5)  in  Eq.(B-2),  the  pressure  becomes 


iNPfU  (r)  ’I  *  dMP’0*>  I - - - 

- d?~  f  f— 56 -  'j'  *  a“(P>Js-uA(?> 

o  p  a 


£e'm<6^o>  J(k-a(rM(IK>ei(*-=>kdkdVf>  (  c-6) 


With  the  trigonometric  t rans format  ions  for  0  and  0 

0  0 


.  .  'NPfUA^  ”  ’’l  af  (p,e  )  - - - —  -irv\o  im(0L cosy  -S°cos9  ) 

PT(x’r’<?o)  ■  —3—  2  I  J  — dQ  °  ■  nI1^  »S(P)P3  »ft(P)e  '""e  »  “b  “ 

2Tr  m=-!=0  o  p  a 

m=2N 


r  p  .  9, 

a  b 


J  Vk-a(r)m)(lK)me 


■  1  /  o  Gi  iL^b  b  .  . 

ka(r)  a(p)J  ( r)  C0Sra  a(p)  COSvyJ 


dkdc  do 

Qf 


(C-7) 


The  ©^-integral  is  as  before  (see  Eq . 1 7  of  text) 
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A^(k-a(p)mje£/a(p) J  =  j 


dfT(p,3a)  i  ( (k-a(p)m)  8P/a ( p) )  cos8 


00 

o  a 


“  dG  (c-8) 

a  '  7 


After  folding  the  m-series  to  m=0  to  +»  ,  the  pressure  of  Eq.(C-7)is 
brought  to 


■NPfU  (r)  _ _ _ _ 

PT(x,r,cpJ  =  - - -  J  Vl  +  a2(p)pP  •  U  (p) 


lx?  p 

"  '  p  ikVif^-  "  ifpy  •  7(77  coscpJ 


I  JA(l<eb/a(p)]k(IK)0e 


r  p  9, 

..a  a  b 


dk 


M  ® 


i(k-a(r)m)(ar-0^coscpQ,)/a(r)  -  i(k-a(p) m)a  p/a(p) 


m5)2N  J('K)nLAV(l<-a(p)m)0P/a(p))(k-a(r)m;e 


+A^(k+a( p)m) 0P/a(p)^k+a( r)m Jt 


(k+a(  r)m)  (or-6^cos9a)  /a  (  r)  -i  (k+a(p)m)of 


/a  ( P)-  \ 

dkj>  dp 


(C-9) 


Let  k-a(p)m=u  in  the  first  term  and  k+a(p)m=u  in  the  second  term  of  the 
second  k-integral.  Finally,  Cq,(C-9)  can  be  written(for  p  <  r)  as 


PT(x,r,cpo)  =  -  —  fV-1  l  U  (p)  V1  +  a2(p)p2 
tT  p 

(  “  r  -  /  .  \  iu(ar-e’'coscp)/a(r)  -  i  uop/a  ( p)  , 

•  <  J  u  1q( up)  Kq( ur)  ^ I  mPa  rt^A^ u0jV a (  p)  Je  e  j  du 

I  o  ~>J 


M  a>.-  .  n 

+  ^  Ji_u-Va(r)-a(p)/jmJ|m(  |u  +  a(p)m|p)K  (|u  +  a(p)m|r) 
m=N,2N  o 


ImPart  ^A^u0P /a  ( p)  Je 


i[u-(a(  r)  -a(p))mj(or-0^coscp  )/a(r)  -iuop/a(p)_, 


du 


(c-10) 

|_Cont 1 d  3 


R-2063 


M  00  - 

+m=?,2N  {  Lu+(a(r)-a(p))lm(|u-a(p)mjp)Km(|u-a(p)ni|r) 

*  {lmPartj^A(u0^/a(p))e'[U+(a(r)  ”a(p))m:i(CTr~6bcosc^a^  /a(r)^-' ucjp/a(p) 
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